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Abstract. We study a family of approximations to Euler's equation depend- 
ing on two parameters e, rj > 0. When e = r) = we have Euler's equation and 
when both are positive we have instances of the class of intcgro-differential 
equations called EPDiff in imaging science. These are all geodesic equations 
on either the full diffeomorphism group Diff^foo (IR* 1 ) or, if e = 0, its volume 
preserving subgroup. They are defined by the right invariant metric induced 
by the norm on vector fields given by 



/ {Le,t)V, v) dx 



where L StV = (I — 2— A) p o(I-^Vo div). All geodesic equations are locally 
well-posed, and the L £ir) -equation admits solutions for all time if n > and 
p > (n + 3) /2. We tie together solutions of all these equations by estimates 
which, however, are only local in time. This approach leads to a new notion 
of momentum which is transported by the flow and serves as a generalization 
of vorticity. We also discuss how delta distribution momenta lead to'Vortex- 
solitons", also called "landmarks" in imaging science, and to new numeric 
approximations to fluids. 



In Arnold's famous 1966 paper [2], he showed that Euler's equation in 1" for in- 
compressible, non- viscous flow was identical to the geodesic equation on the group 
of volume preserving diffeomorphisms for the right invariant L 2 -metric. This raises 
the question, what are the equations for geodesic flow on the full group of diffeo- 
morphisms in various right invariant metrics? Arnold also gave the general recipe 
for writing down these equations but, as far as we know, geodesies of this sort were 
not specifically studied beyond the 1-dimensional case, until Miller and Grenander 
and co-workers introduced them into medical imaging applications. In 1993 they 
laid out a program for comparing individual medical scans with standard human 
body templates [17]. Subsequently they introduced a large class of right-invariant 
metrics on the group of (suitably smooth) diffeomorphisms using norms on vector 
fields given by: 



(Lv, v)dx. 



Here L is a positive definite self-adjoint differential operator. They proposed to 
measure the distance from the subject scan to the template by the length of the 
L-geodesic connecting them (see their survey article [IE])- The geodesic equation 
for these metrics are integro-differential equations called EPDiff (or 'Euler-Arnold' 
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equations). In this paper we want to study the relationship of Euler's equation to 
EPDiff. 

To be specific, we shall use in this paper the group Diff (K n ) of all diffeomor- 
phisms <p of the form ip(x) — x + f(x) with / in the intersection H°° of all Sobolev 
spaces H s , s > 0, and also its normal subgroup Diff,s(]R n ) where / in the space S 
of all rapidly decreasing functions. See |16j for Lie group structures on them. 

Note that, in the H°° case, /, along with its derivatives, will approach as 
||x|| — > oo, but not necessarily at any fixed rate. The geodesic equation in these 
metrics is similar in form to fluid flow equations except that it involves a momentum 
m(x, t) = Lv{x, t), called momentum because it is dual to velocity in the sense that 
the 'energy' can be expressed as \\v\\ 2 L = f(v,m)dx. 

The geodesic equation EPDiff of interest in this paper is this: 

(1) <9 t m = — (v ■ V)m — div(w)m — to ■ (Dv)* 

In coordinates, we can write the right hand side more explicitly as: — ^2j(vjd Xj rrii + 
d Xj Vj ■ rrii + rrijd Xi Vj). Note that v can be recovered from to as v — K * m where 
K is the (matrix-valued) Green's function for the operator L, that is, its inverse in 
the space S. 

The rather complicated expression for the rate of change of momentum - that is 
the force - has a simple meaning. Namely, let the vector field v integrate to a flow 
tp via 

d t ip(x,t) = v((p(x,t),t) 
and describe the momentum by a measure-valued 1-form 

to = rriidxi ® (dx\ A • • • A dx n ) 

i 

so that H^Hf, = / (v, m) makes intrinsic sense. Then it's not hard to check that the 
second equation simply says that to is invariant under the flow ip, that is: 

m(-,t) = (p(;t)*fh(',0), 

whose infinitesimal version is the following, using the Lie derivative (see [T31 3.4]), 

(2) d t m{-,t) = -£„(. )t) m(-,t). 

Because of this invariance, if a geodesic begins with momentum of compact support, 
it will always have compact support; and if it begins with momentum which, along 
with all its derivatives, has 'rapid' decay at infinity, that is it is in 0(||x|| _ ") for 
every n, this too will persist. This comes from the lemma: 

Lemma. [16j If ip £ Diff^-oo (W l ) and T is any smooth tensor on R" with rapid 
decay at infinity, then <p*(T) is again smooth with rapid decay at infinity. 

Moreover this invariance gives us a Lagrangian form of EPDiff: 

(3) dt<p{x,t) = J K*-' t \x,y){<p{y,t)*m{y,Q))dy = K^U (¥>(-,t).m(-,0)) 

where K v (x,y) = K(ip(x),(p(y)) 
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The main result of this paper is that solutions of Euler's equation are limits of 
solutions of equations in the EPDiff class with the operator: 

(4) L e ,„ = (I- ^A) p o (I- ^Vodiv), foranye>0,77>0. 

We will show that all solutions of Euler's equation are limits of solutions of these 
much more regular EPDiff equations and give a bound on their rate of convergence. 
In fact, so long asp > n/2+1, Trouve and Younes have shown [22] that these EPDiff 
equations have a well-posed initial value problem with unique solutions for all time. 
Combining our result with theirs gives a new way of approximating solutions of 
Euler's equation by solutions of a more regular equation. 

Our approach is closely related to several strands of work reformulating Euler's 
equation in a Hamiltonian setting. The first goes back to P. H. Roberts' 1972 paper 
[20j on weakly interacting vortex rings in R 3 where a finite dimensional Hamiltonian 
system for a finite set of such rings is introduced (his equations (35) and (36)). In 
1988 OseledetspjJ] gave a completely general Hamiltonian reformulation of Euler's 
equation. He introduced the dual momentum variables m(x, t) described above, 
called j(x,t) in his paper, that have non-zero divergence in general. With a suit- 
able Hamiltonian, he recovers Euler's equation as a Hamiltonian system. He notes 
that when the momenta are sums of delta functions, one recovers Roberts' system. 
Subsequently, in a second direction, Alexandre Chorin and his students Thomas 
Buttke, Ricardo Cortez and Michael Minion developed the discrete approximation 
by vortex rings as an effective way to solve Euler's equation numerically in line 
with the general "Smoothed Particle Hydrodynamics" (SPH) technique (see [HH]). 
Chorin discusses this technique in §1.4 of his book [5] where he calls the momentum 
variables 'magnetization'. Finally there is a third strand connected to our work. A 
key point in EPDiff is the use of operators L of the form (/ — A) p which have the 
effect of smoothing the velocities v that solve the equation. The case p = 1 arose 
earlier from the study of the Camassa-Holm equation [5], also called the a-Eulcr 
equation for incompressible flows in dimension bigger than one. The CH equation 
is very explicitly related to EPDiff in Holm and Marsden's 2003 paper |9], which 
strongly motivated the present paper. 

The main point here is that all this work, in both the discrete and continu- 
ous cases, fits in logically as special cases of the general EPDiff setup and thus 
as geodesic equations on the group of diffcomorphisms with Riemannian metrics 
depending on two auxiliary parameters. Besides these formal connections we give 
what we believe are new existence theorems for certain cases of EPDiff that, as 
stated above, lead to explicit bounds on the convergence of the particle methods to 
solutions of Euler's equation. In the last section, we show that Roberts' dynamics 
of vortex rings is the same as our geodesic dynamics when the momentum is a sum 
of delta distributions. In this context, it is interesting that this dynamical system 
generates in many cases higher order singularities in the infinite time limit and we 
illustrate these. 

The authors wish to thank Alexandre Chorin, Constantine Dafermos, Darryl 
Holm and Andre Nachbin for very helpful conversations and references. 
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1. OSELEDETS'S FORM OF EULER'S EQUATION 

Oseledets' Hamiltonian formulation of Euler's equation states that for a suitable 
kernel K, Euler's equation becomes equation ^ above. To describe his result, 
consider how EPDiff might be related to Euler. First of all, it's natural to take K 
to be the identity matrix times a delta function 5[x) because then ||u||f is just the 
kinetic energy J \ v(x)\ 2 dx where | • | denotes also the Euclidean norm in K™. Then 
v = 771, and EPDiff becomes: 

d t V = ~[v ■ V)7J - div(7j)7J - V(i|7j| 2 ) 

which looks like Euler's equation if the divergence of v can be made to be zero for 
all time and the last term can be interpreted as the gradient of pressure. But how 
do we keep the divergence of v zero? In fact, the right link between Euler and 
EPDiff is a little more subtle and requires the ansatz: 

Kij(x) = 5ijS(x) + d Xi d x .H 

with the Hessian of an auxiliary function H. With this form of K, we get: 



3 \ j 



or 

v = m + V(ff * div(m)). 
Substituting this into EPDiff and assuming div(u) = we again get Euler's equa- 
tion: 

= 0(777 + [v ■ V)?77 + div(u)777 + 777 • (Dv) 

= d t v - V(dt(H * div(m))) + (v ■ V)u - (v ■ V)(V(ff * div(?77)) + 

+ V(i|w| 2 ) - V(i? * div(m)) • [Dv? 
= d t v + (v ■ V)w + Vp 

with the pressure 

p = -d t [H * div(m)) + \\v\ 2 - [v ■ V)[H * div(m)). 

But now we can also guarantee that the divergence of v is zero if we choose H 
correctly. We have: 

div(Tj) = div(777) + A (if * div(777)) = div ((S Q + AH) * m) 

so all we need to do is to take H to be the Green's function of minus the Laplacian 
and, at least formally, we get Euler's equation. But K now has a rather substantial 
pole at the origin. In fact, if V n is the area of an (77 — l)-sphere, then: 

H(x) _{ (^W(i/N- 2 ) if»>2, 
U[x) \ ^io g (i/M) if n = 2 



so that, as a function 



a a tt 1 nXiXj - Sij\x\ 2 

d Xi d X] H — — ■ — , lf ^0. 

Letting Mq denote this matrix-valued function, note that convolution with any of 
its elements [Mo)ij is still a Calderon-Zygmund singular integral operator defined 
by the limit as e — > of its value outside an e-ball, so it is reasonably well behaved. 
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As a distribution there is another term; compare with |10l Thms. 3.2.4 and 3.3.2]. 
It is not hard to check that: 

a a tt distribution , , , % 1 , , u 

o Xi a X] H = (Mo)ij d {x)6ij 

Now convolution among distributions is associative and commutative so we have 

rn + <9 2 (i/)distr * Tn = rn + VH * div(m) 

which is the identity if div(m) = and has values with div = 0, i.e. is a projection 
onto the subspace of divergence-free vector fields. As it is self-adjoint we see that 

m ^v= (m + <9 2 (iJ)distr) = • m + M * m) 

is the orthogonal projection Pdi v =o of the Hilbert space of vector fields m onto the 
subspace of divergence free vector fields v follows easily. This is the vector form of 
the Hodge decomposition of 1-forms and is orthogonal in each Sobolev space. The 
matrix given by the value of Mq at each point x £ M. n has M.x as an eigenspace with 
eigenvalue n — 1 and Rx 1 - as an eigenspace with eigenvalue —1. So if we let u x be 
the unit vector in direction x and be its orthogonal complement, then we have 
the useful formula: 

Pdiv=o{m){x) = 

2^- • m(x) + • lim / -r-r^ ((n - l)P Uy (m(x - y)) - P u ± (m(x - y)) dy. 

V n e ^J\y\>e \V\ v ' 

With this choice of K, EPDiff in the variables (v, m) becomes the Euler equa- 
tion in v with pressure given by an explicit function of m and v. Letting fh = 
'Y^, i midxihe the 1-form associated to the vector field m, this gives us Oseledets's 
form for Euler 's equation: 

v = Pdiv=o(™) 

(5) 

dtm = — [y ■ V)m — m • (Dv) 

Since v is divergence free we can use fh instead of J^i m idXi (8) dx\ A . . . dx n . In 
integrated form, we have: 

d t <p = Pdi V =o{m) o <p 
(6) ^ ^ 

m(-,t) = <p(-,t)*m(-,0) 

This form of Euler's equation uses the variables v,m instead of the traditional 
v,p (velocity and pressure) but it has the great advantage that m, like vorticity, is 
constant when suitably transported by the flow. In fact, m determines the vorticity, 
defined in arbitrary dimension as the 2-form w = d(J2i v idxi). This is because v 
and m differ by a gradient, sow = dm also. Thus the fact that vorticity is constant 
along flows is a consequence of the same fact for the momentum 1-form fh. This way 
of writing the velocity field as a convolution with a momentum field means we write 
the velocity field as a superposition of the elementary vector fields Pdiv=o(nT,od Xo ) 
for all points xq and vectors m®. In dimension 2, xq — (0,0), mo = (1,0), this is 

(2 _ 2 o \ 
x \ x \i , J with a singularity at where it has a 

dipole as vorticity. In dimension 3, this vector field is an infinitesimal vortex ring 
which is how Roberts' paper [5D] connects to our paper. 
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One of the motivations for this formulation of Euler's equation is that if v(x 1 0) 
are any initial conditions for velocity, we take any momentum m(x,0) such that 
v(-,0) = Pdiv=o( TO ('i 0)). As Chorin has pointed out, in many situations one can 
start with m(-, 0) of compact support and then m will remain of compact support 
even though v will have heavy tails due to the effects of pressure far from the support 
of m. This seems to be one of the reasons why his numerical vortex dipole/vortex 
ring technique works so well. 



However, the above equations ([5]) and ^ are not part of the true EPDiff frame- 
work because the operator K — Pdiv=o is not invertible and there is no correspond- 
ing differential operator L. In fact, v does not determine m as we have rewritten 
Euler's equation using extra non-unique variables m, albeit ones which obey a con- 
servation law so they may be viewed simply as extra parameters. The simplest way 
to perturb this K to make it invertible is to replace the above Green's function H of 
the Laplacian by the Green's function H £ of the positive definite operator e 2 I — A 
for some constant e > (whose dimension is length -1 ). The Green's function may 
be given explicitly using the 'K' Bessel function via the formula 



This has exactly the same highest order pole at the origin as K did and the second 
derivative is again a Caldcron-Zygmund singular integral operator minus the same 
delta function. The main difference is that this kernel has exponential decay at 
infinity, not polynomial decay. By weakening the requirement that the velocity 
be divergence free, the resulting intcgro-diffcrential equation behaves much more 
locally, more like a hyperbolic equation rather than a parabolic one. 

Note that here K e scales as K e (x) — e n Ki(ex) and that, as e goes to zero, 
the limit of K e is just our previous kernel K . Taking the Fourier transform and 
inverting, we can find the corresponding operator L £ in several steps: 



2. Approximating Euler with EPDiff 



H e {x) = c n e n - 2 \ex\ x - n l' 1 K nl2 _ 1 (\ex\) 



for a suitable constant c„ independent of s(see [T]). 
Then we get the modified kernel 




Now the inverse of this as a matrix is the remarkably simple: 




and this comes from the differential operator: 

L e = I — ifV o div 
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Thus we have inverse operators as required by the EPDiff setup: 

v = K e * m, to = L £ (v). 
Finally this operator L E defines the simple metric: 

(y, L e (v)) dxi A • • • A dx n = / (|i>(x)| 2 + 4ydiv(7j)(x) 2 ) dx\ A • • • A dx r , 



As in Arnold's original paper, formally at least, solutions of EPDiff for this 
K e , L £ are geodesies in the group of diffeomorphisms for this metric. EPDiff is the 
geodesic equation with momentum and velocity but in this case it simplifies to a 
form involving only velocity that closely resembles Euler's equation. Substituting 



the formula for L s , we calculate as follows: 



d t (Vi) = {K E )ij * dtirrij) 

= -{K e )ij * {{v k m j>k + rrij ■ div(v) + m k ■ v ktj ) 
= -(K e )ij * ((v k Vj tk + vjdiv(v) + v k ■ Vk,j) 

- ir(^fediv(«),jfe + div(«) i:! div(u) + div(v)^v k ,j)) 
= -{Kjij^vMk + \ {^v{x)f - (^) 2 ) + Vj div(v) - j,(v k div(v),kh 
= -(K e h * (v k v j<k + i (\v(x)\ 2 + (^)) 2 ) . + (L e ) jk (v k div(v))) 
= -(Keia * (v k v j>k ) - Vi dW(v) - \{K e ) l3 * {\v{x)\ 2 + (^l) 2 ) , 

Now we also have the identity: 

(K e * V/), = f,t + d l d J H e * f d = + AH e * f,i = e 2 d Xi H s * f 

3 

so the final geodesic equation is: 

(7) d t (v) = -(K e ) * ((v ■ V)«) - v ■ dW(v) - ^VH e * (\v{x)\ 2 + (i») 2 ) 



This is certainly the simplest choice for a metric which allows non-zero diver- 
gence but, as e — > 0, seeks to make the divergence smaller and smaller so that, in 
the limit, the divergence must be identically zero and we have the L 2 metric on the 
group of volume preserving diffeomorphisms. At the same time, the above equation 
approaches Euler's equation. We will show below that solutions of the above equa- 
tion must approach solutions of Euler's equation and, when the momentum has 
rapid decay at infinity, we will give an explicit bound on the rate of convergence. 
Curiously though, the parameter e can be scaled away. That is, if v(x,t),m(x,t) 
is a solution of EPDiff for the kernel Ki, then v(ex, et), m(ex, et) is a solution of 
EPDiff for K e . 

The above case of EPDiff still has a singular kernel K e for which existence 
theorems are difficult (see below) . The cases of EPDiff which have been analyzed 
and used in medical imaging applications [TH E21 E3] involve kernels which are 
C . We can easily make our singular example a limit of better behaved examples. 
The simplest way is to compose the above operator L e with a scaled version of 
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the standard regularizing kernel (I — A) p giving the positive definite self-adjoint 
differential operator given above (equation (j4j of the Introduction): 



(/ _ i y o div) 



Here the constant 77 has dimension length and although e and ?y could be scaled 
away by themselves, the composite kernel has a dimensionless parameter rj-s. Since 
L e ,r) is a composition, so is its inverse and hence the kernel is now the convolution: 

K e „ = * K £ 

where G^ is the Green's function of (I — ^-A) p and is again given explicitly by 
a 'K'-Bessel function d p ^ n rj~ n \x\ p ~ n ^ 2 K p _ n / 2 {\x\/rj). The reason for inserting p in 
the denominator of the coefficient is that for p ^> 0, the kernel converges to a 
Gaussian with variance depending only on 77, namely (2^fiTr])~ n e~ ][X][ ^ 4v . These 
approximately Gaussian kernels lie in C q if q < p— (n+ l)/2. So long as the kernel 
is in C 1 , it is known that EPDiff has solutions for all time [22] ■ A particularly 
simple case is when p = (n + 3)/2. Then the Green's function is just a constant 
times the C 2 function (1 + \x\/r))e~\ x \/' n as you can verify by taking n = 1 and 
checking that that this is the Green's function of 1 — rj 2 (d/dx) 2 . 



\ \ I / 
^ \ \ \ 



s \ I 
\ \ I 
/ / / 



/ t t t 
I t t \ . 



\ \ \ V 
/ I \ \ 



Figure 1. The dipole given by the kernel Ko tV in dimension 2. 



Finally we may also consider the limiting case e — 0,77 > 0. In this case v — 
*-Pdiv=o ( m ) so v has divergence zero. There is no L because m can be any vector 
field. However EPDiff in Oseledets's form form ([5]) makes perfect sense. Like Euler's 
equation it gives geodesies on the group of volume preserving diffeomorphisms. As 

always, the energy is E = J v ■ m and this is conserved on geodesies. Even though 

2 

we have no L, we can rewrite the energy using (/ — —A) p v — Pdiv=o( TO ) , giving 
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us: 

E = I Vm = J v P div = {m) = I "v(I- ^A) p v. 

The case p = 1 has been introduced and studied by Holm and collaborators (see 
[5], equation (8.29)) who use the letter a for our r\ and call EPDiff the a-Eulcr 
equation: 

(1 - a 2 A){d t v) = -(v ■ V)(l - a 2 A)v - (1 - a 2 A)v ■ {Dv f + Vp, 
div(v) = 0. 

You can also drop incompressibility and when n = 1 this becomes the Camassa- 
Holm equation [5]. 

The K for the e — 0, r/ > metric is just the convolution G ?) * Pdiv=o- This K 
can be explicitly calculated using the fact that the Green's function H is harmonic. 
We use: 

Theorem 1. LetF(x) — f(\x\) be any integrable C 2 radial function on 1" . Assume 
n > 3. Define: 

H F (x) = f min (rj^, j-^) F(y)dy = / F{y)dy + f T^dy 

Jr- Vl ' m > Fl" 2 J\y\<\x\ J\y\>\x\\y\ n 2 

Then Hp is the convolution of F with , is in C 4 and: 

di {H F ){x)=-{n-2)^ [ F(y)dy 

Fl J\y\<\x\ 

d i9j (H F )( X ) = (n-2) ( miX >-J?W 2 f w< m*y ~ V ^ F ^) 

If n = 2, the same holds if you replace l/|a;|"~ 2 by log(l / ) and omit the factors 
(n — 2) in the derivatives. 

Proof. The idea is to first check that Hp £ C 2 with the above expressions for 
its first and second derivatives. This is straightforward when x ^ 0. Near 0, let 
F{x) = a + b\x\ 2 + o(|a;| 2 ). Then one checks that: 



F(y) , \ {n-2)aV n , |2 (n - 2)bV n 



hence the expressions for the first and second derivatives extend across the origin. 
Taking the trace of the matrix of second derivatives, one finds that; AHp = —(n — 
2)V n F. Since the Green's function of —A is l/(n — 2)V n \x\ n ~ 2 , this implies that 
Hp is the convolution F * (l/|a:|"" 2 ). □ 

This applies to F = G, ; for example, or to the limiting case where F = a 
Gaussian, giving the following expression for the kernel -Ko,?j for finite p or the 
limiting Gaussian case: 

(8) K 0>n (x) = 5 jM^l G ^{x) + XiXj ~y |2/n Mean Bw (g^) 
where B a is the ball of radius a centered at the origin. 
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We can summarize all possibilities in a handy table (we have changed notation 
slightly to use double subscripts e,n for all cases): 



no L K ,o = P d iv=o = SijS + didjH 

no L K 0tV — G)j * Pdiv=o _ see above 

L e>0 = I - prV o div K Efi = Sijd + d l djH e 

Ls, v = (l - £ A) P o (7 - £ V o div) K £lV = <%G, ( ; p) + (G^ * H e ) 



3. Existence theorems for the L s>v metric 



It is well known that local solutions of Euler's equation itself, that is Lq^o, exist, 
e.g. see [TTJ [2T]. Moreover global solutions of the EPDiff equations L E ^,e, r\ > 
0,p> (n + 3)/2 have been shown to exist by Trouve and Younes (unpublished but 
apparently implicit in the results of [22] for geodesies in what they call 'metamor- 
phosis'). Their result extends easily to the EPDiff equations Lo^ v because the kernel 
Ko,n is still C 1 , which holds so long as p > (n + 3)/2. The method here is based 
on the Lagrangian form ([3]) of EPDiff. For completeness, we include the proof: 

Theorem 2. Let rj > 0,p > (n + 3)/2 and K — K e ^ be the corresponding kernel. 
For any vector-valued distribution mo whose components are finite signed measures, 
consider the Lagrangian equation for a time varying C 1 -diffeomorphism ip(-,t) with 
ip(x, 0) = x: 

dt<p(x,t) = J K(ip(x,t) - V (y,t)){D V {y,t)y^ T mo(y)dy. 
Here Dtp is the spatial derivative of (p. This equation can be solved for all time t. 



Proof. The Eulerian velocity at (p is: 

V v (x) = J K(x-^(y))(Dip(y))- 1 - T m Q (y)dy 

and W v (x) — V v (tp(x)) is the velocity in 'material' coordinates. Note that because 
of our assumption on m , if <p is C 1 , then V v and W v are C 1 vector fields on K™. 
The equation can be viewed as a the flow equation for the vector field tp ^ W v on 
the union of the open sets U c — {tp <E C 1 (R")™ : dct (Dtp) > c}, where c > 0. We 
claim this vector field is Lipschitz on each U c : 

\\W Vl -W V2 \\ c i <C.\\p 1 -p 2 \\ c i 

where C depends only on c. This is easy to verify using the fact that K is uniformly 
continuous. As a result we can integrate the vector field for short times in C . But 
since (Dtp(y,t))~ 1:T mo(y) is then again a signed finite M™-valued measure, 

J V v( ., t) (x)(Dtp(y,t))- 1 > T m (y)dx = \\V v{ . tt) \\ LsiV 

is actually finite for each t. Using the fact that in EPDiff the L E))7 -energy || V v (.j\ \\ l e ,„ 
of the L 6i?) -geodesic is constant in t, we get a bound on the norm || V v (. t t) de- 
pending of course on r\ but independent of t, hence a bound on || V v t. t) Wc 1 ■ Thus 
||y>(-,i)||co grows at most linearly in t. But d t Dtp = DW V = DV V ■ Dtp which 
shows us that Dtp grows at most exponentially in t, hence det Dtp can shrink at 
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worst exponentially towards zero. Thus for all finite t, the solution ip(-,t) stays in 
a bounded subset of our Banach space and the ODE can continue to be solved. □ 

For L Sf Q, e > we proved in a previous paper j!4j that the L ej o-metric defined a 
well behaved Riemannian metric on the group of diffeomorphisms of K™ in that the 
infimum of path lengths joining two distinct diffeomorphisms was positive. Here 
we prove that for all s and 77, including e = and/or 77 = 0, geodesies exist locally 
- though as in the Euler far as we know, they might become singular in 

finite time hence not be indefinitely prolongable - and that these local solutions 
behave continuously in the parameters e, 77. In particular, as e, 77 — > they approach 
solutions of Euler's equation. 

Everything depends on proving a Sobolev estimate for the time derivative of 
certain energies. We need the following straightforward lemma whose proof we 
omit: 

Lemma. If 77 and e are bounded above, then the norm 

NIL,7,= E / (D a L^v,D a v)dx 

\a\<k 

is bounded uniformly above and below by the metric: 

\\v\\l k + ^\\dW(v)\\l k + J2 r > 2(H ~ fe) / \ D ° V \ 2 + M^ div M| 2 

fc+i<H<fc+ P 

(where H k is the k th order Sobolev norm). 

Assuming k is sufficiently large, for instance k > (n + 2p + 4) works, we now 
prove the main estimate: 



\9 t (IMILJ I < c-IMIlU, 



where, so long e and 77 are bounded above, the constant C is independent of e and 
77. 

Write M V = (I- ^A)p, so that m l = M 11 v l - ^M v div(v) 4 . Using EPDiff and 
integration by parts, the time derivative is given by: 

h 9 t(\Hl, s , v ) = E / (dtD a m,D a v) 

\a\<k JR " 

= E J {—l)^ +1 D 2a Vi.(vj.mij + mi.Vjj+rnj.Vj t i) 

Next replace the m% by M v Vi — -pM v div(v)^. Integrating the third term by parts 
to move the i th derivative of Vj to the other factors and noting that the two terms 
involving the second derivative of divw cancel, one checks that the estimate can be 
reduced to 6 terms all of the form J D 2a f ■ g ■ M v h with one of the triples: 

(f,9>h)=(vi,Vj,Vi,j), (v^v^Vjj), (ui,divu,Ui), 
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Next we expand M v to |/3 1 <p( — l)' /9 ''^pr^ 2 ^ an d integrate by parts some more, 
moving a D a from the first to second or third terms and a D" from the third to 
first or second terms, getting terms 

„2|/3| f 
p\P\ J 

where a = ai + a 2 , (3 — (3i + (3 2 . Now either a.\ or a 2 is less than or equal to k/2 so 
that the corresponding (second or third) term in the integrand has order at most 
k/2 + p + 1, hence < k — n/2 — 1. Thus by the Sobolev inequalities, its sup norm 
is bounded by its k th Sobolev norm and we have: 



p\f>\ 



< 



WnWrP+hfWv, ■ \\g\\ B „ • \\r,^D^h\\ L 2, if \a 2 \ < k/2, 

h m D a +^f\\ L 2 ■ \\ri^D^+^g\\ L 2 ■ \\h\\ Hk , if |ai| < k/2. 

The first term is always bounded by ||/||fc, E ,J7 and so is the other D-term except in 
the first case with a 2 = 0, \a\ = k, \f3\ = p and h is a first derivative with h cither 
a component of v or div v/e. In this last case, the third term has k+p+l derivatives, 
so the lemma does not apply. But we can still integrate by parts, putting the £ th 
derivative on the other terms. If |/?2| > or if / = vg, this reduces again to terms 
bounded by the (fc, e, rj) norm. The only remaining case is when f = h, and then 
we have: 

J D^f ■ g ■ D a +f>U = \ J Vr+Pfft -9 = -\J (D a +fiff ■ g, e 

and this finishes the proof of the estimate. 

Using this estimate, we can prove: 

Theorem 3. Fix k,p, n with p > n/2+l,k > n + 2p+4 and assume (e,r)) 6 [0, M] 2 
for some M > 0. Then there are constants to, C such that for all initial conditions 
vq(x) £ ffk+p+i ^ th ere i s a un iq U e solution v £ . v (x,t) of the above case of EPDiff 
(including the limiting case of Euler's equation) for t G [0, to]. The solution v e ^ v 
depends continuously on e,g £ [0,C] 2 and satisfies \\v £<ri (-,t)\\k,e,n < C f or °^ 
t£[0,t o }. 



We follow a standard approach, used, for example, in [3T], Ch. 16 and 17. First 
consider existence. For e,rj > 0, existence has been proven in [35]. But by our 
estimate and Gronwall's lemma, we have a local upper bound uniformly in e, w for 
these solutions: 

IK„lkfe<c(t),*e [o,to]- 

But, for k,p as above, the Hilbert space with this norm is compactly embedded in 
C 1 (IR™) in the local sense that any bounded sequence for the former has a subse- 
quence which, for every compact subset K C R", converges in C 1 (AT). Therefore 
v e,r)(t) he in a 'locally' compact part of the Banach space of C 1 functions of (x,t). 
Therefore, as e or r\ tend to zero, they have a convergent subsequence whose limit 
v must be a solution of the corresponding EPDiff. Therefore by Gronwall's lemma 
again the original estimates gives H k bounds on this solution. 
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Next we prove that this solution is unique. Let us temporarily abbreviate Lq,e 
by L and let v and v be two solutions of EPDiff for this L. We write u — v — v for 
their difference and follow the ideas of the preceding estimate to estimate 

J^H^Il!, = 2 J (■^Lu,u)dx so using summation of indices: 

= —2 / Ui ■ (vjLvij + Lvi divv + LvjVjj — VjLvij — Lvi divv — LvjVjj) 



= —2 / Ui (ujLvij + VjLui.j + Lui divv + Lvi divu + LvjUj^ + LujVjj) 

Next replace all expressions of the form Luk by Uk — -^dkdivu. Then integrate by 
parts by the "div" part of the last term, that is replace — U{ ■ j^(divu) jVjj by 

Ui • -p(divu).ijVj + divu ■ -p(d\vu)jVj 

The term with the second derivative of divv cancels the term with the second 
derivative of divv arising from the second term VjLuij in the above expression. 
With this and further integration by parts, we get: 

= -2 / UiUj(Lv t j + Vj ti ) + ufdivv + UiUijVj + u l u j . i Lv j + 



u, 



divu Lvi + -ii(divu) y i(divuvi — Uidivi>) 



UiUj(Lvij — Lvj t i + Vj t i) + M^(divu — ^divv)+ 

+ Ui divu (Lvi — Lvi + -p-div(v) 7 i) + (divu) 2 (divv — ^divw) 
< C.\\u\\l 

where the constant depends on the strong sup bounds we have for v and v. By 
Gronwall again, this means that we have a growth estimate on \\u\\ 2 L as a function 
of t. In particular, if u is zero at time 0, it is always zero and this proves uniqueness. 

Finally, as e goes to zero, we again have the solutions lying in a 'locally' compact 
part of C 1 (as above) so if there is only possible limit, they must converge to this 
limit and are continuous in r\. Likewise, as e converges to zero, this solution must 
converge to that of Euler's equation. 



4. Conserved quantities: linear and angular momentum 



We would like to derive the conservation laws from Noether's theorem using the 
fact that our geodesic equation is invariant with respect to the Euclidean group 
SO(n) k R™, as we did in our earlier paper [15] . However, if we take (X, w) € 
so(n) kK™ to be the infinitesimal generator for the 1-parameter group (exp(tX), tw), 
composition maps a diffeomorphism tp G Diff^R") to the diffeomorphsm exp(tX)o 
if + tw. Unfortunately, the latter diffeomorphism no longer rapidly falls towards 
IdjR» so it is not in Diffs(R™). The infinitesimal generator for this action is 

C(x, t0 )(^) = <9t| (exp(fcf) o ip + tw) = X o ip + w. 

For any geodesic t M- <f(-,t) the right invariant inner product G !p (C,^x,w)( l p), ^t) 
should be constant in t, according to Noether's theorem in the form of [31 sec- 
tion 2.6], if the action above was a left action of the motion group on Diffs(R n ). 
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Nonetheless this can be achieved by taking Diffs(M™) as the normal subgroup of 
an extension of the motion group which can be described as a group of diffeomor- 
phisms which fall rapidly to "Euclidean motions near infinity". Instead of doing 
this in detail we directly check that the the above well defined expression is indeed 
constant in t along each geodesic. Note first that 

G<p{C(x,w)(<P),<Pt) = G V (X o ip, (fi t ) + G v (w, (ft 

= G U (X, ipt o tp- 1 ) + G u (w, (p t o tp- 1 ) 

= G u (X,v) + Gi d (w,v) where <p t {-,t) = v(<p(-,i),t) 



= J(X.x,L(v)(x))dx + j \w,L(v)(x))dx 
= J(X.x, m{xj) + J (w, m(x)); 



the first expression viewed as a linear functional in X € so(n) is the so(n)*-valued 
angular momentum mapping. If we identify so(n)* with so(n) via the Killing form 
we can write the angular momentum succinctly as J xAm(x). Similarly the second 
expression leads to the linear momentum given by J rh(x). 

Let us finally prove that these momenta are conserved by the geodesic flow. We 
shall use the geodesic equation in the form dtfh = —L v fh. Then we have 

d t J (X.x,fh(x)) = J(X.x, d t m{x)) = - J(X, C v m) 

= J(C v X,m)= J([v,X],m) = J(-C x v,m) 

= — J (£xv,Lv)dx now use Cx{L) = and Cx{dx) = 0, 

= - J \C x {{v,Lv)dx) = 0. 
For the linear momentum the proof is similar. 



5. Explicit bounds on the approximation I 



Assume you start with the same initial condition v(x, 0) and integrate with both 
Euler's equation and EPDiff with L e ^. Exactly how close are they? If you look 
at the kernels K e>ri , you see that the effect of e > is to shrink the tails of K 
from polynomial to exponential and correspondingly, to eliminate the pole of its 
Fourier transform at zero. On the other hand, the effect of 77 > is to smooth the 
singularity of K at zero or to suppress the high frequencies in its Fourier transform. 
These being opposite operations, we need to estimate their effects separately. In 
this section, we consider the case 77 = and compare Euler's equation with that 
given by L £) o- Let vo(x,i) be the solution of Euler's equation and let v e (x,t) be 
the solution of EPDiff with L e , (below abbreviated to L £ ). Our goal is to prove 
the theorem: 

Theorem 4. Take any k and M and any smooth initial velocity v(-,0). Then there 
is are constants t ,C such that Euler's equation and {e 1 Q)-EPDiff have solutions 
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vq and v e respectively for t € [0, to] and all e < M and these satisfy: 

IM-,t)-tfe(- s t)|| ff * < Ce. 



Note that by Theorem [3] we have essentially any bound we need on both Vq 
and v e . The to is needed only to guarantee the bounds on the solutions derived in 
Theorem [3] hold for a big enough k to give us the needed bounds. As above the 
proof is based on an estimate of the form: 

(9) d t (\\v - v £ \\ 2 Hk ) KCi.Wvo-VeWlt+eCiWvo-VeWu* 

where G\ and C2 depend on the initial condition v(-,0) and k but not on e at all 
times t for which Theorem [3] holds for needed norm bounds on vo and v E . 

Let u = vq — v e and calculate as follows, using the geodesic equation (J7J) for L e : 
hdt(\\uf Hk ) = J2 I D a u-{D a {d t vo-d t v e )) 

\a\<k'' 

= £ / D a u ■ D a { - Ko * (v ■ V)v 

|a|<fc 

+ K e * (v e • V)w £ + v e • div(w £ ) + ^(VH e ) * \v e {x)\ 2 + |(Vff e ) * (divi; e ) 2 } 

< ]T PNU- {|p a (i^ e -X )*(«e-V) V£ || L2 + p Q (z; e -div(i; £ )|| i2 

|a|<fe 

+ \\D a (^(VH s ) * \v E \ 2 )\\ L2 + ||D«(|(VF £ ) * (divt, e ) 2 )|| L2 

+ \\K * D a {{u ■ V)« e )|U>} - £ f D a u-K Q * D a {{vo • V)u) 

|a|<fc 

< ||u||^ k • {||(# * (« e • V)u e || ff * + |K-div(« e )|| H * 
+ ||^VH £ * \v £ \ 2 \\ Hk + U 2 VH e * (div Ue ) 2 ||^ 

(10) +||((^_ t , 6 ).vK|U»}- J] f D a (v -K *v e )*D a ((vo-V)u) 

\a\<k 

Here, in the last line, we used the fact that i^o, being an orthogonal projection, 
has norm 1 and is self-adjoint. Likewise K e has a Fourier transform that at fre- 
quency £, is a matrix with eigenvalues 1 and e 2 /(e 2 + |£| 2 ), hence is also a bounded 
self-adjoint operator with norm 1. 

For the first term, if we abbreviate v £ to v, first write: 
((K - K e ) * (v ■ V)u) = F - (Ko - K e ) * v ■ div(u) 

where F t = 2j(ifo - K s ) tJ * (vjVk),k = 2J 0fc(i£o - ^e)y * WjVfc 

3,k j,k 

The Fourier transform of the derivative of the difference of the K's is: 
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Thus 



VjVk\ 



i j,k 

< /(i+ia fc kP (4iw Efe^) 2 -Z 

< /a + iei 2 ) fc ^E i« 2 = tE n^«*iii* 

j,k j,k 

Thus wc find: 

\\(K - K e ) * {v ■ V)v)\\ m < | E ll«i"fc||ff* + ||«- div(»)|| fffc ) < Ce\\v\\ 2 k+m>efi 

for some universal constant C coming from the product rule for Sobolev spaces and 
m = fn/4] . 



The second term in the expression ()10|) has a similar bound. To treat the third 
and fourth terms we need the bound on the norm of convolution with \7H £ : 



hence ||eVff e * f\\n k < ill/llir* f° r anv function / and in particular: 



\(%rVH e ) * \Vs(x)\ 2 \\ Hk < mVe\ 2 \U < Ce\\v £ \ 2 



2 » "ej • rer; uh* — 2 1 

\2|| . s' e II / diyjis 



|(iVi/ e )*(div U£ ) 2 || ff , < 



Iff* < Cte||«. 112 



k+m,e,0 

ellfe+m,e,0 



For the fifth term in expression (|10j) we want sup bounds on k + 1 derivatives of 
u e and so we can get a bound: 

||(« • V)v e \\ Hk < C\\u\\ h h ■ \\v e \\ H t, with i = k + 1 + rail] 
We come to the last term in (|T0l) . Write it as: 
(11) I D a u- (D l} v Q -V)D a - f} u+ f D a (v e -K *v e )-D a {v -V)u 

0</3<a,\a\<k 

In the first term of dill) , the summand with /3 = vanishes because it equals 
J wo • V( ) and vq has zero divergence. Using a sup norm on D^Vq, the 
remaining summands are bounded by \\vq — v E \\ 2 Hk times this sup norm. This sup 
norm is bounded by a universal constant times \\vo\\e with I — k + f^^i]- To 
bound the second term in (|lip , using the expression for Kq we find v £ — K * v £ = 
ViJ * div('We). Now calculate: 

J D a (v E - K * v e ) ■ D a ((v • V)«) =-J2j D a (didjH * div(« e )) • £> a ((« );%) 

» J 

But didjHo has Fourier transform (£i£j)/|£| 2 , a matrix with eigenvalues and 1, 
so the L 2 norm of the first factor is bounded by e||w E ||fc je ,o- Then, as above, we get 
a bound of the form: 

Ce||w e ||fc, s ,o ■ \\v \\h> ■ IMIff* 
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with I = k + 1"^^"]. Now using Theorem [3J we see that we can bound all needed 
norms of vq and v £ on this time interval by norms of the initial condition v(-,0). 
Putting everything together, we get the asserted bound. 

To complete the proof of the Theorem, simply apply Gronwall's lemma to dt \\u\\ jjk 
(without the square) and deduce that the difference vq — v £ is 0(e) as required. 

In comparing Euler's equation with EPDiff for (e, 0), a key point is that Kq = 
Pdiv=o and K £ = K £ $ have identical singularities at the origin, but their difference 
is much better behaved. In fact convolution with Kq — K £ equals 

V o (convolution with J £ ) o div 

where J £ has Fourier transform l/|£| 2 (e 2 + |£| 2 ). Near the origin, this looks like 
e~' x ' in R 3 , has a log pole in K 4 and is like l/|x|™ -4 in higher spaces. Considering 
Euler's equation and EPDiff for (e, 0) in Lagrangian form ([3]), they differ only by 
changing the convolution on the right hand side by this term. This makes it seems 
reasonable to conjecture that if solutions of (e, 0)-EPDiff do not blow up, i.e. exist 
for all time, then neither do the solutions to Euler's equation. Or conversely, if 
Euler's solutions do blow up, so do solutions of this EPDiff. 



6. Explicit bounds on the approximation II 

Now we want to compare solutions of EPDiff for e > 0, r\ — with solutions 
for e > 0,i] > 0. The difference here is a convolution with the Gaussian G v , so 
solutions with ij > are essentially just smoothed or low-pass version of those with 
7] = 0. We will prove: 

Theorem 5. Take any k and M and any smooth initial velocity v(-, 0). Then there 
is are constants t ,C such that {e,0)-EPDiff and (e,rj)-EPDiff have solutions Vq 
and v v respectively for t € [0, tg] and all e,rj < M and these satisfy: 

ING,*) - «j7(/)*)IIh* < Crj 2 . 



A basic tool is the simple estimate: 
(12) \\f-G^*f\\ L2 <^\\Af\\ L2 

To prove this, just take Fourier transforms and use the elementary inequality: 

(i-(i+^it p ) <v 2 \e- 

Working as in the setup of the Theorem, let too and m n be the momenta corre- 
sponding to vq and Vrf. Write u — vq — v n and calculate the time derivative of: 

HL,o= E [(D a u,D^L e u) = \\u\\ 2 Hk + \\^\\l k . 

H<fc 



18 



DAVID MUMFORD AND PETER W. MICHOR 



We get a lot of terms: 

\a\<k 

(13) = Y ! D ° u ■ ° a {( y o ■ V ) m o ~ K ' V ) m ^ + ( J - G n } ) * K ■ V)(m n ) 

+ div(v )m — div(v v )m v + (I - G&>) * (div(v t) )m,) 
+m ■ (DvoY ~ m v ■ (Dv n )* + (I- G%>) * m v • (£>«„)*} 

By the bound (|12p. the three terms with / — Gjf are bounded by ||w||_fffc times 

2 2 2 

2*pr||0vV)m ); || J j*+a and ^s||m, • div(u,)|| H *+a and ^r||m^ • )*|| Hfc +2. Hence 
if £ = 1 + max(fc + 2, [■§]), then, by the product rule for Sobolev norms, all three 
terms are bounded by Cr] 2 ■ • ||m^||^ for some constant C depending only on 

k and n. Using Theorem [31 this is bounded by C'r] 2 , where C is another constant 
now depending on the initial data as well as k and n. 

If u = toq — Triq, we can write the remaining terms in (|13[) as: 

(w-V)m 7) , (vq-V)u, m r) div(u), udiv(^), m v -(Du) t , u-(Dvq) 

Next use the calculation: 

u = m - G^ * m v -(I- G ( n p) ) * to„ 

= L e (vq — Vn) + term bnded by 77 2 A(m r) ) 

= u — p-V(div(u)) + term bnded by 77 2 A(to 7) ) 

The 77 A (m^) terms are bounded like the previous ones. We finish the proof by 
applying the same tricks we have seen above to the remaining terms. Letting C 
denote suitable constants depending on bounds for vq and v^, the terms with u, 
not u, are easy: 

J2 [ D a u ■ D a ((u ■ V)m„) < C\\u\\ 2 Hk 

Y J D a u- D a (dWu-m ri ) < C\\u\\ H k\\divu\\ H k < C\\u\\ 2 keS) 

Y J D a u-D a {m v -{Du) t ) = -J2 J D a div(u)-D a {u-m v ) + D a u-D a {u-{Dm ri ) t ) 

a a 

< C||«||^* + C||«|| H *||divtt|| H * < C\\u\\l £fi 

Finally, the u terms have two more pieces, one where it is replaced by u and the 
other with -^Vdiv(u). If it is replaced by u, everything is bounded as above by 
Cljull 2 ^,. but where the usual trick is needed: 

J D a u ■ D a {{v ■ V)tt) = J D a u ■ (v ■ V)D a u) + terms with VD p u, (i < a, 

the latter being bounded by and the former being equal to 



{v -V)\D a lu\ 2 = -\ \ div(v )-\D a u\ 2 . 



ON EULER'S EQUATION AND 'EPDIFF' 



19 



The div terms have the 1/e 2 factor but also a cancellation and reduce to: 

- i J div^o)!^"^! 2 + terms with V^^,/3 < a. 



and 



- / D a Ui ■ D c 



a ( div(u^) div(u),i 



D Q div(u) • D° 



a ( div(u^) div(u) 



D 



a UiD a ^ div ("'7).' div(u) \ 



which have the needed bounds. Finally, using Gronwall's lemma, the Theorem is 
proven. 

7. Approximating Euler solutions via landmark theory 



The great advantage of having a C 1 kernel is that we can now consider solutions 
in which the momentum m is supported in a finite set {Pi, • • • , Pn}- The support 
is called the set of landmark points and in this case, EPDiff reduces to a set of 
Hamiltonian ODE's based on the kernel K = K e , v , e > 0, r) > 0: 



Energy E = ^ m al Kij(P a - P b )m jt 



a . b 



dP a 

dt 

dm a . 
dt 



= Y^K lJ {P a -P b ) 



m bj 



m a jm bk 

b,j,k 



where a, b enumerate the points and i,j, k the dimensions in MP. These are essen- 
tially Roberts' equations from [20]. His paper takes n = 3 so that the landmark 
points are the center of 'circular vortex rings'. He assumes they do not get too close 
to each other and takes K(x) at all x 7^ to be the Euler kernel Pdiv=o 5 our -K"o,o- 
He sets K(0) = SijK for a constant k which comes out of the specific model used 
for each finite (non-infinitesimal) vortex ring. What using our kernel Kq.jj does is 
just smoothly interpolate between the kernel Pdiv=o at points x far from - but 
which is singular at - and a C 1 function near with K(0) — Sij.K. 

For some other PDEs (like the KdV or Camassa-Holm equations) solutions whose 
momenta are sums of delta distribution are called solitons. In analogy to this we 
can call vortex- solitons or vortons the solutions with momenta supported in finite 
sets. 

All of our kernels have the form Kij(x) = G\(\x\)b~ij + d Xi d Xj G2(\x\) hence at 
every point x have eigenspaces Ka; and (Rx) . For any vector x, let x = p x ■ u x , 
where p x — \x\ and u x is a unit vector; and let P Ux be the projection to the subspace 
R • u x and P u ± be the projection onto the perpendicular subspace (R • Ux)- 1 . Then 
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the matrix Kij (x) can be written in terms of two scalar functions K\ and K 2 as 

K(x) = K x {p x )P u „ + K 2 { Px )P ui , if x ± 

and as n ■ I at the origin. If K\ = K 2 , then Kij would be a multiple of the identity 
and we would have the case studied in our previous paper |12) . But this never 
happens for our metrics. For example, in the Kq^ case, using formula ([8]) and the 
fact that Grf 1 is a monotone decreasing function of \x\, we get: 

Ki(x) = |Mean B|x| (G^)(x) 

K 2 (x) = Gf(x) - lMean Bj jG^)(x) 

(Jfi - K 2 )(x) = Mean Bb| ((l7W)(x) - G%\x) > 0, 

k = |G^(0) > K x {x) > K 2 {x) if x ^ 

If we differentiate the formula for K, we get the following formula for its deriv- 
ative: 

D v K(x)=K'- L (\x\){v,u x )P u „+K 2 (\x\){v,u x }P u ± 

K x {\x\) -K 2 (\x\), ± . , ± ± , 

\x\ ~^( V ' U ^ ® U x + u t ® U x ) . 

Using this we can rewrite the geodesic equations in a geometric form: 

E = l mQ | 2 + K i(Pab)(Pu ab m a , P Uab m b ) + K 2 (p ab )(P u ± b m a , P u ± b m b ) 



dPg 

dt 

dm a 
dt 



a a^b 

b m b + K 2 (p ab )P u ± m b 

by&a 

(K' 1 (p ab ){P Uab m a ,P Uab m b ) + K 2 (p ab )(P u ± b m a ,P u ± b m b )j u ab 



b^a 

K l(Pab) ~ K 2 {p ab ) / \ 

> \{m a ,u ab )P u ±m b + {m bl u ab )P u ±m a 



a#6 



One of the characteristics of these landmark space EPDiff geodesies as that when 
two landmarks near each other, they can either repel or attract. If their energy is 
low compared to their angular momentum, they repel and vice versa. When they 
attract, they typically spiral in towards each other with the momentum of each 
landmark point growing infinitely while their sum remains bounded. They do not 
collide in finite time. Whether this characteristic reflects developing singularity 
behavior in Euler's equation is not clear because, as soon as landmarks approach 
closer than rj, solutions of EPDiff are no longer close to those of Euler. This 
attraction is clear with only two landmark points but, at least in the case of the 
Weil-Peterson metric on cosets of Diff(5' 1 ), following a geodesic typically produces a 
hierarchical clustering of many landmarks (unpublished work of Sergey Kushnarev 
and Matt Feizsli). 

We want to look at the simplest cases of one or two landmark points. One 
landmark point is very simple: its momentum must be constant hence so is its 
velocity. Therefore it moves uniformly in a straight line I from — oo to +oo. As a 
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geodesic in Diff (R n ), it will push everything in front of it, compressing points ahead 
of it on t while pushing out points near I to maintain incompressibility. Behind the 
landmark, they will be sucked back towards i to compensate for the rarihcation left 
by its passage. By rotational symmetry around t and time-reversal symmetry, the 
motion, from t = —oo to t = +00 c an only be a shear in which points are dragged 
forward parallel to i by a distance which goes to zero as you go further from I and 
goes to 00 as you approach £. 




Figure 2. The result of the incompressible flow from t = —00 
to t = with n = 2, momentum concentrated at one point and 
K = Ko,r,, V =■ 



Now consider the case of two landmark points Pi , P2 with momenta mi , rri2 ■ 
By conservation of total momentum, mi + m,2 is a constant. We can reduce this 
Hamiltonian system by fixing the total momentum m and dividing by translations. 
We get a new system in the variables SP = P2 — Pi and 5m = rri2 — mi with 
equations of motion: 



E = k ■ (jSm\ 2 + |?ti| 2 ) + Kij(Sp) ■ (vriiVfij — SniiSmj) 



1 .1 



djSPt) 
dt 

d{8rrii) 
dt 



■ Srrii ~ Kij(SP) ■ 5m j 
3 

} J d Xl K ]k {8P) ■ (mjlrTk - SmjSmk) 
jk 
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or, letting 5P = p ■ u for a unit vector u, in geometric form: 

E = k- {\5m\ 2 + \m\ 2 ) + K^p) {\P u m\ 2 ~ \P u 8m\ 2 ) 

+ K 2 (p) (\P u ±m\ 2 - \P u ^Sm\ 2 ) 

d(6P) c . . „ r 

= k ■ 5m — Ki(p) ■ P u 5m — K 2 {p) ■ P u ±5m 



dt 

"' ( " m) = - (K[(p)(\P u m\ 2 - \P u Sm\ 2 ) + K 2 (p)(\P^m\ 2 - \P^5m\ 2 ) u 



dt 



— 2 -2— {{m, u)P u ±m ~ (<5m, u)P u ±5m) 



Note that the derivatives of Sp and 5m lie in the span of 5P, 5m and m. Thus 
this three dimensional space is constant in time so we can assume SP, 5m, m € M 3 . 
The total angular momentum is: 

uj = Pi A mi + P 2 A m 2 = Pl + P2 A m + \5P A 5m. 

If m = 0, then the two vectors 5p, 5m always lie in a fixed two dimensional space 
and their cross product is constant, equal to 2w. We can then make a further sym- 
plectic reduction and compute what happens in terms of the three scalar variables 
p, (5P,5m), \5m\ which moreover must lie on one sheet of a hyperboloid: 

4|w| 2 + (6P,5m) 2 = p 2 ■ \5m\ 2 , p ■ \5m\ > 2|w|. 

The energy then simplifies to 



E = (k — K 2 (p))\5m\ 2 - w-W (6P, 5 



in : 



"-(5P,5m} 2 +^ 



UJ 



Its level curves on the hyperboloid must then be the geodesies. Note that as long 
as the kernel is C 2 , (n — Ki(p))/p 2 and (Ki(p) — K 2 (p))/p 2 are finite at the origin 
hence bounded. 

We can illustrate this in the simple case of 3-space with kernel Kq,i,P — 3. As 
stated above, then the smoothing kernel is C 2 and has the elementary expression 

(x) = (1 + \x\)e-^ = l-Ml + ... . 

It's easy to calculate the mean of this function over a ball and we get: 

Mean B|x| (g^) (x) = 24|x|- 3 (l - e"M (l + |z| + ^ + ^)) 

OO \ 

n=4 / 




hence 

K^l-\x 2 + --- , K 2 = l~lx 2 + --- . 
A typical plot of the contours of E in the (p, |<5m|)-plane is shown in Figure |31 



In the figure, if an orbit hits the heavy black line defined by p ■ \5m\ = lo, 
then (5P,5m) is instantaneously zero and, along its orbit, changes sign. On the 
two-sheeted cover given by including this sign, this is a smooth orbit in which p 
decreases to a minimum where (5P, 5m) = and then increases. One sees that 



ON EULER'S EQUATION AND 'EPDIFF' 



2:S 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

IS m| 



Figure 3. Level sets of energy for the collision of two vortons with 
m = 0,?7=l,cj = l. The coordinates are p and |<5m| and the 
state space is the double cover of the area above and right of the 
heavy black line, the two sheets being distinguished by the sign of 
(8m, 5 P). The heavy black line which is the curve p ■ \5m\ = uj 
where (6m, SP) = 0. Each level set is a geodesic. If they hit 
the black line, they flip to the other sheet and retrace their path. 
Otherwise p goes to zero at one end of the geodesic. 



there are two types of orbits: scattering orbits where the vortons separate infinitely 
at both t — ±oo and p has a minimum at some point in time; and capturing orbits 
which either start or end at infinity but spiral indefinitely, getting closer and closer, 
at the other limit. Which happens depends on the relative size of the angular 
momentum and the energy exactly as in the simpler case studied in |12j . Here if 
E > (8/5)|w| 2 , the points attract while if E < (8/5)|o;| 2 , they scatter. 

When the landmark points attract, this simple system forms higher order sin- 
gularities. If we take coordinates so that 8P is on the cci-axis and 5m in the 
(x\, X2)-plane, then for p very small, we have: 

P 1 = (p/2,0,0), P 2 = (-,9/2,0,0), mi = -(0,0), m 2 = —(0,0) 

P P 

where 2C = (SP, 5m), hence (using the limiting values of the fc-terms in the formula 
for energy) we get C 2 « \E— 2uj 2 > 0. Then, as these points approach each other, 
the corresponding global vector field in M 3 approaches: 

(K Oil ) l .(x + (p/2,0,0))-(K Otl ) l .(x-(p/2,0,0)) (C 

Vt[X) ~ p '{ LJ 

This vector field for uo = 2C is illustrated in Figure [H Whereas for any column 
vector A, Kqi ■ A is a vortex ring with maximum norm at the origin and maximum 
vorticity along a ring centered at the origin and lying in a plane perpendicular to A, 
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its derivative v is now zero at the origin and it has maximum vorticity there. In our 
case, computing the derivatives Dv(0), we find that near the origin, the flowlines 
of v spiral in along the (xi,X2)-plane and shoot out along the a^-axis. 




Figure 4. Streamlines and MatLab's 'coneplot' to visualize the 
vector field given by the xi-derivative of the kernel i^o.i times the 
vector (1,2,0). See text. 



Another case which is easy to explore is when m lies in the plane spanned by SP 
and 5m. The angular momentum no longer descends to a function on the SP, 8m 
space but we may numerically integrate the geodesic equations. Figure [5] shows 
geodesies all starting with the same SP and 8m but with varying m fixed along the 
y-axis. 

It is extremely easy to compute landmark geodesies numerically even in much 
more complex situations and we hope that, letting 77 — >■ 0, this may be a useful tool 
to exploring the instabilities of Euler's equation itself. 
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